    autoPtr<radiationModel> radiation(radiationModel::New(thermo.T()));
